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Abstract 



It is well known that sedimentary rocks having same porosity can have 
very different pore size distribution. The pore distribution determines many 
characteristics of the rock among which, its transport property is often the 
most useful. Multifractal analysis is a powerful tool that is increasingly used 
to characterize the pore space. In this study we have done multifractal anal- 
ysis of pore distribution on sedimentary rocks simulated using the Relaxed 
Bidisperse Ballistic Model (RBBDM). The RBBDM can generate a 3 - D 
structure of sedimentary rocks of variable porosity by tuning the fraction p 
of particles of two different sizes. We have also done multifractal analysis on 
two samples of real sedimentary rock to compare with the simulation studies. 
One sample, an oolitic limestone is of high porosity ( 40%)while the other 
is a reefal carbonate of low porosity around 7%. 2 — D sections of X-ray 
micro-tomographs of the real rocks were stacked sequentially to reconstruct 
the real rock specimens. Both samples show a multifractal character, but 
we show that RBBDM gives a very realistic representation of a typical high 
porosity sedimentary rock. 

Keywords: multifractal, sedimentary rocks, simulation, pore size distribu- 
tion 
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1 Introduction 

Sedimentary rocks are often the storehouses of natural oil and gases whose 
extraction depend on the permeability of these fluids through them. The 
transport properties of sedimentary rocks depend not only on the porosity 
of the rocks but more importantly on the pore size distribution (PSD) and 
their connectivity. The pore space can be a continuum of of pores with ex- 
tremely varying pore sizes ranging over a scale of 10^, besides being extremely 
complex and heterogeneous and often self-similar. 

Fractal and multifractal analysis are increasingly used to study complex 
heterogeneous systems which show self-similarity on several length scales. 
They have the ability to provide an accurate representation of the hetero- 
geneous pore geometry and address the relationship between porosity and a 
range of physical processes happening in a porous medium like transport of 
water in soils, extraction of oil and natural gases and CO2 sequestration in 
sedimentary rock. Multifractal analysis has been done using fractal models 
(Rieu and Sposito, 1991), image analysis of two-dimensional sections of soil 
blocks (Tarquis et al.,2003; Dathe et al.,2006; Grau et al.,2006), analysis of 
three-dimensional pore systems reconstructed by computer tomography (Tar- 
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quis et al.,2007), mercury intrusion porosimetry (Vidal Vazquez et al.,2008) 
and nitrogen absorption isotherms (Paz Ferreiro et al., 2009). 

In this work, the authors use the Relaxed Bidisperse Balhstic Deposition 
(RBBDM) to simulate a three dimensional porous rock structure of varying 
porosity and pore distribution. In our efforts (Giri et al., 2012a, 2012b) 
to probe the geometry of the microstructure of the pore clusters produced 
by the RBBDM at different porosities, the authors had noticed that the 
simulated structure had a fractal nature over different length scales. The 
power law exponent had different values over different length scales which 
hinted that the pore space might have a multifractal nature. This was further 
strengthened by diffusion studies in connected pore clusters. For the entire 
range of porosities studied, diffusion was found to be anamolous with different 
values of diffusion exponent over different length scales. Real rock samples 
were studied for comparison of simulation results, and similar signature of 
multifractal nature was found there! We shall therefore investigate whether 
our simulated porous structure generated with the RBBDM is indeed a 
multifractal in its pore distribution. We shall bring out the differences in 
the PSD of the simulated structure at different porosities through a study of 
their multifractal spectral dimensions. Finally we shall compare our results 
with similar studies done on real limestone and carbonate rock samples. 

The details of RBBDM have been discussed in earlier works by the au- 
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thors (Sadhukhan et al. 2007a, 2007b, 2008, 2012) in the study of various 
transport properties like permeability and conductivity through sedimentary 
rocks. A brief outline of the model will be given here for the sake of com- 
pleteness. The basic algorithm is to deposit particles of two different sizes 
ballistically. In 3-D (2 + 1 model), we drop square 1x1x1 and elongated 
2x1x1 grains on a square substrate. It is well known that natural sand 
grains are angular and elongated (Pettijohn, 1984) , so the aspect ratio 2 
is realistic. The cubic grains are chosen with a probability p and elongated 
grains with probability (1 —p). The presence of the longer grains leads to 
gaps in the structure. The porosity (p, defined as the vacant fraction of the 
total volume, depends on the value of p. For p = 1, a compact structure is 
produced. As p is decreased, isolated pore 'clusters' start appearing and the 
porosity increases. For a specific value of p, the threshold value, a structure 
spanning cluster is generated. However this cannot be called a 'percolation 
threshold' as in the case of random percolation. In this respect the RBBDM 
is different from random percolation problem (Stauffer and Aharony, 1994). 
The RBBDM being a modification of the Random Deposition Model, for 
p = 1, the surface width keeps increasing with height. In fact in the limit 
of infinite height, at least one narrow structure spanning pore cluster is al- 
ways present. When p is gradually decreased to below 1, larger grains are 
introduced and the grains settle on the structure following the Ballistic De- 
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position Model. The introduction of even an infinitesimal quantity of the 
larger grain, can close a deep surface trench creating an elongated pore clus- 
ter. Obviously these pore clusters are longer near the surface of the structure 
than at the bottom. The presence of a single large grain sitting atop a long 
pore cluster, introduces correlation between adjacent columns (Karmakar et 
al.,2005). As the fraction of large grains increase, the correlation spreads 
through the system. So a substrate of sufficient height needs to be generated 
before the porosity value can stabilize. Our model is different from the ran- 
dom percolation problem as even an infinitesimally small fraction of larger 
grains introduces correlation between columns, thus robbing the system of 
its randomness. 

The RBBDM has the potential of generating a structure with a con- 
nected rock phase that is needed for any stable structure, and a tunable 
porosity. As the fraction of longer grains is increased, unstable overhangs 
can develop. If a larger particle settles on a smaller particle, a one-step over- 
hang is created. If a second larger particle settles midway on the previous 
large particle, a two-step overhang is created if there is no supporting parti- 
cle immediately below the protrusion of the second overhang. This two-step 
overhang is not stable and the second large particle topples over if possible, 
according to the rule scheme as shown in fig.([T]). This leads to to compaction. 
In their earlier works (Manna et al.2002; Dutta and Tarafdar, 2003), the au- 
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thors have shown that the sample attains a constant porosity only after a 
sufficient number of grains (depending on sample size) have been deposited to 
overcome substrate effects. Here, a Lx x Ly x Lz size sample was generated, 
from which a Lx x Lx x Lx sample was selected after the porosity had sta- 
bilized to within 0.001 percent. The selected sample was chosen from below 
the deepest trough at the surface to eliminate surface effects. All simulation 
was carried out on this sample. To check for finite size effects, we carried out 
our studies for Lx = 32,64,128,256 for which Lz = 1000,2000,4000,7000 
respectively. The results reported in this work did not show any finite size 
dependence. All results on simulation are reported for 256 x 256 x 256. 

Figs. ([2^) and ([2]d) show vertical sections, (x-z) and (y-z) planes of the 
generated sample at maximum porosity (pmax = 0.42. Fig. (lit) shows a hor- 
izontal section (x-y plane) of the sample at the same porosity value. The 
anisotropy in the pore geometry is clearly visible. The pore clusters have an 
elongated and interconnected appearance along the z-direction while the dis- 
tribution of pores along the horizontal plane is quite homogeneous. As the 
fraction of larger grains is decreased, the porosity of the sample decreases 
and the pore distribution becomes more anisotropic nature. Fig.([3K) and 
fig-dSb) show the vertical, (x-z) and horizontal sections, (x-y)), respectively 
of the sample at = 0.073, a very low porosity. The elongated isolated pore 
clusters are prominent in the direction of assembly of the grains, whereas the 
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pores remain homogeneously distributed in the (x-y)plane. 

To compare our simulation results with real rock samples, X-ray tomog- 
raphy micrographs of 2 — D sections of two real sedimentary rock samples 
obtained from an oolitic limestone (pure calcite) from the Mondeville for- 
mation of Middle Jurassic age (Paris Basin, France), and a reefal carbonate 
from the Majorca Island, Spain, have been used. The oolitic limestone is 
composed of recrystallized ooliths with a mean diameter of less than a few 
hundred fim. Each pixel of both the micrographs corresponds to 5.06 micron. 
For every real rock sample studied, each micrograph section was converted 
to a binary file form such that corresponded to a pore site and 1 corre- 
sponded to a rock site. The binary file was then converted to a grey-scaled 
picture, as shown in fig.(j3]), using MATLAB. An array of 1000 consecutive 
sections were put together precisely to reconstruct the binary file form of the 
real three dimensional rock structures. In each of the two samples, the real 
structure chosen was 1000 x 1000 x 1000 in size, and all study on real rock 
was carried out on this structure. The 2 — D sections of oolitic limestone cut 
in the direction of assembly (growth) from the reconstructed 3 — D structure 
are shown in Figs.(j4^),(x-z)plane, and (Ilb),(y-z) plane. Comparison between 
figs, (jl^) and ([Dd) shows the pore distribution to be isotropic. Fig. (^) shows 
a 2 — section of the bedding plane of the same rock structure. There seems 
to be slight anisotropy in the pore distribution in the bedding plane and the 

8 



direction of growth. Similar 2 — D sections of the reefal carbonate rock sam- 
ple have been cut along the direction of assembly of the reconstructed 3 — D 
sample and shown in figs.dHi and|4^). Fig.(j4]F) shows a.2 — D section of the 
same sample cut along the bedding plane. Unlike the previous sample, any 
anisotropy that may be present in the pore size distribution along the two 
directions, is not easily discernible. 

In the following section we shall briefly summarize the multifractal con- 
cepts and estimation techniques. We shall follow up with our analysis and 
results on the RBBDM. This is followed by our studies on real rock struc- 
tures where we compare our findings with the results of simulation. Finally 
we conclude with a discussion on the importance of such analysis to study 
pore distribution in sedimentary rocks. 

2 Multifractal Concepts 

Highly inhomogeneous systems which do not obey self-similar scaling law 
with a single exponent, may actually consist of several intertwined fractal 
sets with a spectrum of fractal dimensions. These systems are said to be 
multifractal. Such systems have a complex distribution which arises from 
peculiarities of their generation. These are not a simple collection of fractal 
systems, rather one may say that these constitute a distribution of several 
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fractal subsets on a geometrical support. Each of these subsets is character- 
ized by its singularity strength and fractal dimension. 

Usually multifractality arises from a spatial distribution of points, each 
having a strength or weight associated with it. The weights may also have a 
non-trivial distribution. A growing diffusion limited aggregate (DLA) with 
weights proportional to the growth probability assigned to each site, is an 
example of such a multifractal. However, a system with equal weights as- 
signed to each point, may form a geometrical multifractal. The distribution 
of pore clusters in a sedimentary rock belongs to this class. 

The multifractal system has local fractal dimension a{x), which may be 
determined by applying the 'sandbox' method at different points x on the 
system. Different regions distributed over the system may have the same 
local dimension a{x). If we collect the xs with the same a{x)s, we shall have 
identified one fractal subset of the multifractal. The scaling exponent for 
this subset has a value, say f{a). A plot of f{a) versus a gives a typical 
multifractal spectrum. For a monofractal, a would be the same everywhere 
and the /(a) versus a curve would reduce to a point. 

Practically it is however not convenient to determine local fractal dimen- 
sions a{x) and therefore f{a). What is done in practice is determine the q^^ 
moments of distribution of points on the system, with q varying from —oc 
to +00. The mass exponent Tq corresponding to the measure of moment q 
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has a one to one correspondence to a{x) through its derivative, and through 
a Legendre transformation to /(a). The details of this correspondence is 
discussed in the following section. 



3 Multifractal Exponents 

Multifractal analysis involves the estimation of three functions: mass expo- 
nent Tq, singularity strength (or local scaling index) ag, and multifractal or 
singularity spectrum f{a). 

Multifractal analysis requires that the chosen system of size be divided 
into a set of different boxes of equal size e. A common choice is to consider 
dyadic scaling down, i.e., successive partitions of the system in k stages 
(fc = 1,2,3....) that generate a number of cells A^(e) = 2^ of characteristic 
length e = L X 2~^. The boxes N{e) cover the entire system and each such 
box is labelled i. 

The probability mass function /ij(e) describing the portion of the measure 
contained in the i*'^ box of size e is given by 

rrit 

where rrii is the number of pore sites in the i^^ box and is the total number 
of pore sites in the entire system. A pore site refers to a pixel that is vacant. 
The measure of the q^^ moment in the i*'* box of size e is termed /i^. /ij(e) is 
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the total number of pore pixels inside the i box of size e. Here q can vary 
from +00 to —00. 

The partition function e) for different moments q is estimated from 
yU^(e) values as 

Nii) 

X(g,e) = 5^/x?(e) (2) 

i=l 

The parameter q describes the moment of the measure. The box size e may 
be considered as a filter so that by changing e one may explore the sample 
at different scales. So the partition function e) contains information at 
different scales and moments. The sum in the numerator is dominated by 
the highest value of fii for g > and the lowest value of /ij for g < 0. 

The measure of the g*'^ moment of the mass distribution of the system is 
defined as 

N 

M{q,e) = J2f^1^' = N{q,e)e' (3) 

i=l 

where 

JV 

N{q,e) = J2f^l-^~^^'^ (4) 

i=l 

N 

If ^^/U^(e) in the limit e — > crosses over from to 00 as c? changes from 
a value less than r(g) to a value greater than r(g), then the measure has a 
mass exponent 

d = r{q) (5) 

The measure is characterized by a whole sequence of exponents r(g) that 
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controls how the moments of probabihty /ij scale with e. For multifractally 
distributed measures, the partition function e) scales with e as 

x(g,e)oce--(^) (6) 

The probability mass function e) also scales with e as 

^^^{q,e) = e^^ (7) 

where is the Holder exponent or 'singularity exponent' or 'crowding in- 
dex' of fi peculiar to each i^'^ box. Greater the value of Holder exponent, the 
smaller is the concentration, and vice versa. Singularity exponents of multi- 
fractal distributions show a great variability within an interval {amax, ttmm) when 
e tends to zero. For a monofractal, this interval reduces to a point. 

Again, the number N^{a) of boxes of size e that have a Holder exponent 
between a and a + 6a obeys a power law as 

N,{a) oc e~-^(") (8) 

where /(a) is a scaling exponent of the boxes with a common a, called 
the singularity exponent. A plot of f{a) versus a is called the singularity 
spectrum, /(a) is the fractal dimension of the set of points that have the 
same singularity exponent a. There can be several such interwoven fractal 
sets of points each with its particular value of /(a). Within each such set, 
the measure shows a particular scaling described by a. 
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Following (Chhabra et al., 1989), the functions a and /(a) can be deter- 
mined by Legendre transformation as 

Thus the singularity exponent defined by eq.([7]) becomes a decreasing func- 
tion of q. Larger values of g (g >> 1) correspond to smaller exponents 
and therefore higher concentration of measure. Similarly smaller q values 
correspond to higher exponents and lower concentration of measure. As q 
varies, points define a parabolic curve that attains a maximum 

value /(ao) the point a^. is the mean value of the singularity expo- 
nents and f{ao) gives the fractal dimension of the support as obtained by 
the box-counting method. 

Another equivalent description of the multifractal system is obtained from 
Dq — q plot, where Dg, called the generalised dimension, corresponds to the 
scaling exponent for the g*'^ moment of the measure. It is defined by 

"^-0 ~ 1-q log{e) 
For the particular case of g = 1, eq.f ITU]) becomes indeterminate, and is 
estimated by I'Hopital's rule. Dq is related to the mass exponent r(g) by 

Tiq) = il-q)Dq (11) 

The generalised dimensions Dq for g = 0, g = 1 and g = 2 are known as the 
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Capacity, the Information (Shannon entropy) and Correlation Dimensions re- 
spectively. Mathematically, the multifractals can be completely determined 
only by the entire multifractal spectrum. However a few characteristic func- 
tions may be used to describe the main properties of multifractals. 

4 Multifractal analysis of sedimentary rocks 
4.1 Simulated rock structure 

In the case of the simulated structure, a 256 x 256 x 256 cube was selected from 
below the deepest trough from the surface after the porosity had stabilized in 
an initial structure of 256 x 256 x 3000, and after the porosity had stabilized 
for a particular choice of p. This system was then covered with hypercubes 
of size e = 2^' with k ranging between 1 to 16. The partition function e), 
calculated according to eq.([2]) for different values of box size e and for different 
moment values q was determined. A log-log plot of xiQi versus e when 
plotted, showed a deviation from linearity beyond a certain range of e. A 
power law scaling was observed only in the range e = 2^ to 2^. All calculations 
have been done within this range of e for q ranging from —9 to +9. The 
exponent r(g) for each such q and for every p studied was noted. 

The scaling properties observed in the partition function can be be char- 
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acterized by determining if the scaling is simple as in monofractal, or multiple 
as in multifractal. Figs. ([5^ andlSb) show the variation of r(g) versus q for a 
low porosity = 0.07 corresponding to p = 0.9, and a high porosity = 0.42 
corresponding to p = 0.5 respectively. The data points from simulation stud- 
ies is shown as open circles in both the graphs. The plots for all the other p 
values studied, he within the limits set by these two plots. It is clear that the 
r(g) functions which would have been straight lines for monofractals, deviate 
from linear behaviour. Moreover the slopes of the r(g) for g < are quite dif- 
ferent from those for g > 0. This clearly indicates multiple scaling behaviour, 
i.e. the low density and high density regions of pores scale differently. 

Subsequently, the generalised dimension Dq were estimated in the range of 
q values from -|-9 to —9. Figs.([n^ andinb) show plots of Dq versus q for a low 
porosity corresponding io p = 0.9 and a high porosity value corresponding to 
p = 0.5 respectively. The data points from the simulation study are shown as 
open circles in the figure. For monofractals, all the DqS would lie on the same 
horizontal line. In the case of the simulated sedimentary rock, the first three 
generalised dimensions are different for every p (hence porosity), as shown in 
fig-dZI)- This indicates that the structure is multifractal. Tablel shows the 
results of the generalised dimensions for the first three moments calculated 
for the simulated structure at different porosities. It is apparent that at lower 
p values, between to 0.6, the first three generalised dimensions have values 
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very close to each other. The variation of porosity with p is also very low. 
This indicates that the structure is quite homogeneous here. 

The capacity dimension, Dq provides information about how abundantly 
the measure, defined by eq.©, is distributed over the scales of interest. Ex- 
cept for an abrupt increase at p = 0.7, the value of Dq remains almost the 
same indicating that the same pore abundance is present at all the length 
scales studied. Dq shows a maximum (fig|7j)at p = 0.7. In an earlier study 
(Sadhukhan et al., 2007b) on conductivity through connected pore space 
of sedimentary rocks, the effective conductivity of the simulated rock using 
RBBDM, showed a maximum for p = 0.7 despite a maximum porosity at 
p = 0.5. For structures generated by using the RBBDM , the maximum 
backbone mass of the connected cluster corresponded to p = 0.7. The au- 
thors had established that it was the backbone mass of the connected cluster 
that was most effective for transport. The maximum value of Dq at p = 0.7 
indicates that the capacity dimension is directly related to the pore distribu- 
tion in the backbone of the connected cluster. 

From Tablel, it is seen that the entropy or information dimension Di 
decreases monotonically with increasing p. Lower Di indicates greater con- 
centration of pores over a small size domain, i.e. greater clustering. When 
Di is close to 0, it will be be reflected as a sharp peak on a pore size dis- 
tribution curve. Comparison of figs. ([2] and |3]) clearly indicate that there is 
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greater clustering of pores with decreasing porosity. For the particular case 
of p = 0.7 where conductivity through simulated structures using RBBDM 
showed a maximum, it may be noted that the difference — Di is a min- 
imum. Here the Capacity Dimension Dq is a maximum while the Entropy 
Dimension Di is a minimum, and this has optimized connectivity. This is 
manifested in transport property having maximum values here. 

The correlation function D2 describes the uniformity of the measure (here 
pore cluster size) among different intervals. Smaller D2 values indicate long- 
range dependence, whereas higher values indicate domination of short range 
dependence. From Tablel we see that D2 shows a slight decrease at lower 
porosities which is indicative of long range correlations appearing between 
pores. This is a manifestation of our growth algorithm for the sedimen- 
tary rocks. When the fraction of larger grains is small, elongated and iso- 
lated pore clusters are more prominent. Thus the pore sites show greater 
auto-correlation along these clusters. When porosity increases, though D2 
increases somewhat, it is not too significant as the pore clusters still retain 
their elongated appearance in spite of greater connectivity between the pores. 

The a and /(a) values of the singularity spectrum were computed with 
the help of eq.Q. The plot of /(a) versus a is shown in fig.® for the different 
values of p studied. The shape and symmetry parameters from the singularity 
curves is listed in Tablel for different values of p. The Holder exponent ao 
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for each p gives the average values of local mass distribution for a given scale. 
A greater value of ao indicates a lower degree of mass concentration. This in 
turn would indicate that the pore distribution is highly heterogeneous and 
anisotropic with fluctuations in local densities. From Tablel it appears that 
ao remains almost constant for different porosities of the simulated structures 
except for p = 0.7, where it is a maximum. 

Width of the /(«) spectra is defined as the difference between the a values 
of the most negative moment g_ i.e. amax, and the most positive moment 
g_i_ i.e. amin- The wider the spectrum, i.e. greater the difference between 
(o^max — C(min), the higher is the heterogeneity in the scaling indices of pore 
mass and vice versa. The largest f{a) obtained for g = corresponds to 
the capacity dimension Dq. Small /(«) values indicate rare events (extreme 
values of the PSD). Asymmetry in the /(«) spectra indicate the dominance 
of higher or lower values of pore masses. If the width on the left, — is 
larger, this indicates the domination of large values in the PSD. From Tablel 
we see that ao — Q^mm is larger than ao — ctmax for p values between 0.1 to 
0.6. This implies that there is greater dominance of pore masses here. On 
the other hand, a large right width ao — amax would indicate the dominance 
of extremely small values in the PSD. For the porosities studied, it is clear 
from Tablel that this occurs for p values between 0.7 to 0.9. In this region 
clustering of pore sites into elongated isolated channels leaving larger sections 
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of structure pore free. Fig. ([2]) illustrates this arrangement of pore clusters. 
The /(a) spectra is almost symmetric about /(ao)- This indicates that the 
rock structure has the most isotropic pore distribution here. 

4.2 Real rock structure 

To compare our simulation results with real rock samples, X-ray tomography 
micrographs oi 2 — d sections of two real sedimentary rock samples obtained 
from an oolitic limestone (pure calcite) from the Mondeville formation of 
Middle Jurassic age (Paris Basin, France), and a reefal carbonate from the 
Majorca Islands, Spain, have been used. The limestone is composed of re- 
crystallized oolite with a mean diameter of less than a few hundred /xm. 
Each pixel of the micrographs corresponds to 5.06 micron. Each section 
was converted to a binary file form such that corresponded to a pore site 
and 1 corresponded to a rock site. The binary file was then converted to a 
grey-scaled picture, as shown in fig.(jl]), using MATLAB. An array of 1000 
consecutive sections were put together precisely to reconstruct the binary 
file form of the real three dimensional rock structure. The real structure 
was 1000 X 1000 x 1000 in size. Figs.(|l^ andSb) show sections of the three 
dimensional limestone rock cut along the direction of assembly while fig. (jit) 
shows a section of the same cut perpendicular to the direction of assembly. 
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Similar sections were cut from the reefal carbonate and these are shown in 
figs.(j4ji, andHf). The porosity of the oohtic hmestone was determined 
from the reconstructed rock and found to be 0.073 while the reefal carbonate 
was found to have a porosity of 0.399. This high contrast in their porosity 
values is evident from the panels of fig.Q. 

To compare the results of the RBBDM for sedimentary rocks with real 
rocks, we have plotted the variation of r(g) and Dq versus q for the real 
rocks along with their closest matching porosity samples generated by the 
RBBDM. The log- log plot of r(g) versus q for q ranging from —9 to +9 
for the low porous limestone is shown in fig.(l5^) while the same plot for the 
high porous carbonate sample is shown in fig.([5]D). The non- linear nature of 
the plots with two distinct slopes for positive and negative q values clearly 
indicate that the real rock samples are also multifractal. It is clear that 
at high porosities, the real and the simulated rock give a very good match. 
The simulated rock at high porosities, fig.(|2]), look more isotropic and start 
resembling real samples. For very low porosity like 0.073, the RBBDM 
does not yield a realistic rock sample. The long narrow pore clusters, fig.dS^ 
andlSb), an artefact of the generation rule, are responsible for a pronounced 
anisotropy which results in this mismatch. 

The plots of Dq versus q for the same samples over the same range of q 
values are shown in figs. ([6^ and|6jD) along with their corresponding matches 
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from the RBBDM structures. Once again the similarity between the real 
and simulated rocks at high porosity values is clear. Even for high porosity, 
fig-dSb), the Dq values at higher positive q values show a mismatch between 
real and simulation study. An examination of figs.([2]and|lll) reveal that pore 
cluster size and shape distribution is quite different even though the porosity 
values match. The generalised dimensions corresponding to the first three 
moments in each of the two real rocks, are enlisted in Tahlell and Tablelll. 
It is clear that Dq, Di and D2 are very different from each other in the case 
of the limestone sample showing clear multifractal nature. The difference be- 
tween the first three moments of the carbonate rock though less pronounced, 
is finite. Though both the real rocks have a Capacity Dimension Dq of almost 
the same value, the limestone has a lower value of Di in comparison to the 
carbonate. This indicates greater clustering of pores in the limestone sample. 
One can expect that the limestone will be more efficient for fluid transport 
than the carbonate sample. The smaller D2 value of the limestone indicates 
that there is greater long range correlation between the pore clusters here 
than in the case of the carbonate sample. 

The f{ot) spectra of both the real rocks along with their corresponding 
simulated rock structures having similar porosity, are shown in fig. ([9]). The 
real and simulated rock show similar f{a) spectra at high porosity. At very 
low porosity, the RBBDM fails to create realistic sedimentary rocks. Both 
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the real samples have a wider width of their /(a) spectra indicating that 
there is a greater heterogeneity in the scaling indices of their pore mass. 
With the RBBDM, this nature is observed as porosity increases. It is clear 
from the fig.dHfe') as also from Tahlell that «o — «maa; is greater than Oq — Omm 
in the limestone sample. This indicates that there is greater dominance of 
extremely small values in the PSD. Not only is the porosity of the oolitic 
limestone small, the pore clusters are small and sparse. This is also observed 
in figs.(jlK, andHt). Tahlell I indicates that in the carbonate rock, the 
difference between — amax and — amin is quite pronounced. The larger 
value of ao — amax shows a greater dominance of smaller values in the PSD. 
This dominance of smaller pore clusters is also seen from figs.dHi, andHf). 
Large pore clusters are far and in between here. 

5 Conclusions 

Multifractal analysis on sedimentary rock structures simulated by using RBBDM 
was done. The structures at different porosities, all showed multifractal char- 
acteristics. The complex heterogeneity of the pore size distribution has been 
quantified by the multifractal parameters. The Capacity Dimension gives a 
measure of the pore distribution in the backbone of the connected cluster. 
Fluid transport through such rocks maybe related to the multifractal param- 
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eters Do and Di. A combination of higher Dq and lower Di will result in 
more efficient transport properties. 

Multifractal analysis performed on real sedimentary rock samples showed 
that these too, were multifractal in nature. A comparison of the multifractal 
characters of both the real rocks studied, and their corresponding simulated 
structures with almost matching porosities, was done. The RBBDM showed 
a very good match with the real sample at high porosity values for all the 
multifractal parameters. At very low porosities however, even though both 
the simulated and real samples showed multifractal nature, the character 
match was not so good. At very low porosities the anisotropic nature of the 
simulated structure becomes more pronounced. A more suitable toppling rule 
may perhaps reduce this anisotropy somewhat. In the case of real rocks with 
low porosity, the pore distribution remains more isotropic. We may conclude 
that the RBBDM may be considered a good model for sedimentary rock 
generation especially at high porosity values as it shows similar geometric 
features as real rocks. 
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Table 1: Multifractal parameters for simulated rock. 



D 


(t) 


dn 


di 


d2 












0.0000 


0.4260 


3.0239 


3.0102 


2.9820 


3.0239 


2.8040 


3.2370 


0.2199 


0.2131 


0.1000 


0.4342 


3.0218 


3.0098 


3.0160 


3.0218 


2.8859 


3.2609 


0.1359 


0.2391 


0.2000 


0.4428 


3.0500 


3.0094 


2.9800 


3.0500 


2.8020 


3.3900 


0.2480 


0.3400 


0.3000 


0.4498 


3.0158 


3.0089 


3.0110 


3.0158 


2.8649 


3.3029 


0.1509 


0.2871 


0.4000 


0.4541 


3.0239 


3.0083 


3.0136 


3.0239 


2.8630 


3.3239 


0.1609 


0.3000 


0.5000 


0.4557 


3.0128 


3.0077 


2.9930 


3.0128 


2.8280 


3.3840 


0.1848 


0.3712 


0.6000 


0.4525 


3.0128 


3.0068 


3.0041 


3.0128 


2.8430 


3.5679 


0.1698 


0.5551 


0.7000 


0.4410 


3.1720 


3.0055 


2.9840 


3.1720 


2.8070 


3.4059 


0.3650 


0.2339 


0.8000 


0.4162 


3.0216 


3.0030 


3.0167 


3.0216 


2.6959 


3.4069 


0.3257 


0.3853 


0.9000 


0.3518 


2.9990 


2.9968 


2.9370 


2.9990 


2.6979 


3.7000 


0.3011 


0.7010 
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Table 2: Multifractal parameters for oolitic limestone (pure calcite) from the 
Mondeville formation of Middle Jurassic age (Paris Basin, France). A system 
size 512'^ was coarse grained to 64'^. 






do 


di 


d2 


«o 










0.077 


2.9727 


2.216 


2.6942 


2.9727 


2.0889 


5.1529 


0.8838 


2.1802 



Table 3: Multifractal parameters for reefal carbonate from Majorca island 
(spain). A system size 1000'^ was evaluated. 






do 


di 


d2 


ao 






^0 (-^min 




0.399 


2.989 


2.970 


2.9603 


2.989 


2.898 


3.329 


0.091 


0.34 
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Figure 1: Toppling rule of the larger grains - when a larger grain develops 
a two-step overhang, marked 1 in the figure, with at least two vacant sites 
immediately below the overhang, it topples over in the direction indicated 
by the arrow to assume a more stable state, marked 2. 
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Figure 2: (a) shows x — z section of simulated structure for (pmax = 0.45, 
i.e. high porosity, (h) x — y section at same porosity. Structure looks more 
isotropic. The white indicate pore clusters. 
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Figure 3: (a) and (b) show x — z and y — z section of simulated structure 
for (p = 0.072, i.e. low porosity, matching the porosity of the real rock, z- 
axis indicates the vertical direction. The white indicate pore clusters, (c) 
X — y section at same porosity. Anisotropy in pore cluster structure is quite 
pronounced. 
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XY XZ YZ 




(a) M (c) 




(d) Ce) (f) 

Figure 4: Each section is a square of side 2.58 x 10~^m.(a) and (b) show 
sections of the real hmestone cut in mutually perpendicular planes. These 
are perpendicular to bedding planes, (c) shows a section of the bedding 
plane. These are sections of oolitic limestone (pure calcite)from the Mondev- 
ille formation of Middle Jurassic age (Paris Basin, France), (e) and (f) show 
sections of reefal carbonate obtained from Majorca Islands, Spain, perpen- 
dicular to bedding planes, (f ) shows a section of the carbonate rock along 

the bedding plane. This rock structure looks more isotropic. 
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Figure 5: (a)Variation of r(g) versus q for low porosity (p ~ 0.07. The 
open circles show data for simulated structure while X indicate data of real 
limestone sample, (b) Variation of r(g) versus q for high porosity ~ 0.4. 
The open circles show data for simulated structure while X indicate data of 
real carbonate sample. 
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Figure 6: (a)Plot of D — q versus q for low porosity corresponding to 0.07. 
The open circles show data for simulated structure while X indicate data of 
real limestone sample. (h)D — q versus q for high porosity (f) ~ 0.4. The 
open circles show data for simulated structure while X indicate data of real 



carbonate sample. 
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Figure 7: Variation of the first three generahsed dimension with porosity for 



simulated rock structure. 
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Figure 8: Plot of /(a) versus a is for different porosities of simulated rock. 
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Figure 9: (a)/(a) versus a for low porosity corresponding to </> ~ 0.07. 
The open circles show data for simulated structure while X indicate data 
of real limestone sample. (b)/(a) versus a for high porosity (p ~ 0.4. The 
open circles show data for simulated structure while X indicate data of real 
carbonate sample. 
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